Magnon-mediated thermal transport in antiferromagnets: the link to 
momentum-resolved magnon lifetime data 
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Transport currents in solids decay through coUisions of quasiparticles with each other and with 
defects or boundaries. Since information about coUisional lifetimes is difficult to obtain, most 
calculations of transport properties rely on parameters that are not known independently. Here, 
we use magnon lifetime data for the two-dimensional antiferromagnet Rb2MnF4 to calculate the 
magnon-mediated thermal conductivity without any adjustable parameters, thereby quantifying the 
influence of scattering from domain boundaries on transport. Related strategies have the potential 
to enhance our understanding of thermal transport by electronic and phononic quasiparticles greatly. 
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Thermal transport phenomena are of great importance 
for many areas of science and technology, ranging from 
correlated electron systems jlj 12j to the design of nanos- 
tructured materials for solar energy harvestingiS . A pre- 
requisite for progress towards nanoscale and quantum de- 
vices is thorough understanding of the effects of thermal 
noise and quantum coherence and their relation to size, 
dimension, and coupling to the environment. In quantum 
materials, thermal effects are of fundamental importance, 
as they affect the transport, storage, and processing of 
energy and information. 

It is widely assumed that the thermal conductivity 
mediated by phonons in simple insulators and undoped 
semiconductors is well understood. In fact, even in the 
most advanced investigations of the thermal conductiv- 
ity of elemental semiconductors reported to date, multi- 
ple adjustable parameters describing the collisional rates 
of longitudinal and transverse phonons had to be in- 
troduced to fit the experimental data[lj. Recent ex- 
periments on two-dimensional (2D) antiferromagnetic 
insulators [5HS] have allowed the determination of the 
thermal conductivity Kmag mediated by magnons. Their 
theoretical description requires fewer parameters than 
that of phonons because only two degenerate (or nearly- 
degenerate) magnon branches are present. However, the 
resulting data were analyzed only at low temperatures, 
where Kmag is limited by defects, and the defect density 
could not be determined independently. Clearly, a more 
complete and quantitative understanding of thermal con- 
duction in such model systems is required as a basis for 
theories of complex materials that are of current scientific 
and technological interest. 

We present an approach to this problem that estab- 
lishes a quantitative link over a broad range of temper- 
atures between the macroscopic property of the ther- 
mal transport and the microscopic processes that un- 
derlie it, without the need for adjustable parameters. 



This approach is made possible by the advent of a new 
neutron spin-echo technique by means of which excita- 
tion linewidths can be measured with high resolution in 
energy and wavevector over the full Brillouin zone [9], 
combined essentially with a theoretical description of 
these linewidths that is quantitatively accurate. Very 
recently, magnon linewidths in 2D and 3D antiferro- 
magnets were measured over wide ranges of tempera- 
ture using this technique, and a comprehensive theoret- 
ical description that includes numerical evaluation of 4- 
magnon (two-in/two-out) scattering processes, combined 
with a small contribution from boundary scattering, was 
developedpU]. 

In general, experimental identification of 
ier in magnetic systems of low dimensionality. Since 
the phononic contribution to the thermal conductivity 
K is roughly isotropic, subtracting k measured in a direc- 
tion of negligible magnetic dispersion from that measured 
along a dispersive direction yields an estimate of Kmag 
associated with the dispersive magnons. In a number of 
quasi-lD and 2D magnetic systems, Kmag appears as a 
broad peak that is well separated in temperature from 
the phononic peak[5l [TlJ [12]. As the spinon excitations 
in quasi-lD compounds occupy a continuum in energy, 
measurement of their linewidths is inherently more diffi- 
cult, and corresponding theoretical descriptions are less 
well developed. We will therefore consider only higher- 
dimensional magnetic systems in our analysis. 

The Boltzmann transport equation describes the irre- 
versible approach of a perturbed dilute gas to equilib- 
rium. It can be used to quantify heat flow in the pres- 
ence of a temperature gradient. In order for a perturba- 
tive interaction in a quantum system to allow the system 
to approach equilibrium in the manner described by this 
equation, a large phase space of scattered states must 
exist, and the coupling between scattered states must be 
weak enough that there is essentially negligible correla- 
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tion between them[T3]. The scattering events are then 
elastic, with probabihties given by Fermi's Golden Rule. 
Magnon-magnon and phonon-phonon interactions at low 
temperatures in higher-dimensional systems satisfy the 
above conditions [T3]. Exponential relaxation of such in- 
teracting quasiparticles to equilibrium yields Lorentzian 
excitation lineshapes. A comprehensive description of 
these relaxation rates permits quantification of the heat 
current and the drift distribution of the quasiparticles. 

The heat flow per unit area Q = kVT due to 
thermally-excited collective modes with thermal conduc- 
tivity K in the presence of a thermal gradient VT can be 
calculated from the deviations duk of the mode occupa- 
tion numbers from equilibrium by means of the relation 
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(1) 



where is the excitation energy and its group veloc- 
ity, and Vo is the volume of the unit cell. Values of Suk 
can in turn be calculated from the excitation lifetimes 
by means of the approximate solution of the Boltzmann 
transport equation proposed by Callaway [Mj. Callaway 
considered phonons with linear spin- wave dispersions. In 
the following, we proceed as in his calculation, general- 
izing it to apply to arbitrary spin-wave dispersions. We 
treat magnons, rather than phonons; the factor of 2 in 
the above expression for k originates from the aforemen- 
tioned twofold degeneracy of the spin waves. 

According to the Boltzmann equation, the mode occu- 
pation numbers evolve in response to both 1) the heat 
flowing into a volume element over a local thermal gra- 
dient and 2) changes in mode occupation numbers due 
to collisions. In steady-state heat flow, the mode occu- 
pation numbers remain constant, and thus the change in 
the mode distribution due to collisions is balanced by the 
variation in mode occupation over the thermal gradient. 
This is expressed as 



dnk duk 



(2) 



where Uk is a mode occupation number and duk/dt its 
rate of change due to collisions. Approximation of by 
the Bose distribution Uk (equivalent to the assumption 
that deviations from equilibrium are small yields 



dnk 
dt 



hujk 



(3) 



Callaway was the first to distinguish between the dif- 
ferent effects of normal (i.e. momentum-conserving) and 
Umklapp (i.e. non-momentum-conserving) scattering 
processes on the thermal conductivity. He employed the 
characterization that the equilibrium state approached 
through normal processes reflects the displaced Bose- 
Einstein distribution 1161 



n(A) 
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where 't{T) = fi{T)VT/T lies in the direction of the 
thermal gradient; P(T) is determined below. In contrast, 
Umklapp processes cause the mode occupation numbers 
to relax to the static equilibrium distribution nk- Using 
the relaxation-time approximation |15j, which postulates 
exponential relaxation to equilibrium [IB], Callaway ex- 
pressed duk/dt as 



duk _ n{\) - Hk _^ Uk - Uk 
dt r„ T„, 



(5) 



where r„ represents excitation lifetimes corresponding to 
normal processes and r„ those corresponding to Umklapp 
processes. In Callaway's treatment, these are not simply 
added, as under heat flow the thermal distribution of 
the gas relaxes differently for the two types of scattering 
processes. 

Since A • k is small, n(\) can be approximated to first 
order as 



n{\) « Uk + 



Ak 



(^pht^t/kBT _ ij 



(6) 



The deviations Suk = nk — Uk can then be calculated 
from (5), using (3) and (6): 
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dnk = -TcCk H z — k I • VJ 



hLOk/kBT 



TnflUJk 



where T-1 = r"! r"! Ull 



(7) 



The magnitude of A(T), and thus of j3{T), is deter- 
mined by applying the condition that the normal colli- 
sional processes conserve momentum, namely that 



k 



n{\) - Uk 



0, 



(8) 



which through application of (5), (3), and (7) yields the 
expression 



^ Y.k{rc/Tn){Ck-^T)Ck\^/{hUk) 

Efc(VT«)(l-^c/T„)(k-VT)(7fck/(;iC.fc)2' 



(9) 



where 



Ck 



ksT'^ (^(,hu,k/kBT _ i~)2 



is the heat capacity of a given spin-wave mode. Insertion 
of (7) and (9) into (1) yields n. 

The thermal conductivity calculated using the 
well-known single-mode relaxation time (SMRT) 
approximation[T6] corresponds to /3 = 0, as in this 
approach normal processes are assumed to behave like 
Umklapp ones and relax to the (undisplaced) Bose 
distribution. In Callaway's approach, normal processes, 
since they relax to a displaced equilibrium distribution. 
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FIG. 1. The magnon thermal conductivity of Rb2MnF4 
(black line), computed as described in the text. The grow- 
ing magnitude up to ~ 15 K of the heat current deriving, in 
Callaway's formulation, from the assumption of a displaced 
equilibrium distribution (dashed red line), relative to that 
corresponding to the SMRT approximation (dashed blue line) 
demonstrates the increasing importance of transport by the 
drift current with increasing temperature. Inset: the mean 
free path in the interacting magnon gas (solid line) and the 
mean distance i/2 (dashed line) traveled by a magnon be- 
fore colliding with an intracrystalline boundary with average 
spacing L = 0.58 ± 0.04 umfTO]. 



in fact contribute to thermal drift: a heat current flows 
even in the absence of a temperature gradient |16j. 

In order to use (7), (9), and (1) to calculate the heat 
current, the lifetimes (inverse linewidths) r„ and t„ cor- 
responding to each relevant scattering process must be 
known. Previous applications of Callaway's theory in 
calculations of the thermal conductivity relied on the use 
of adjustable parameters in expressions for r„ and for 
each coUisional process thought to occur [H [T3], as these 
were not known independently. 

For certain types of microscopic coUisional processes, 
including scattering between phonons and scattering be- 
tween magnons, the respective contributions of r„ and Tu 
to the combined lifetime cannot be distinguished by 
means of experiment. However, in a numerical evaluation 
of the linewidths corresponding to a given process, it is 
straightforward to identify the contribution of each com- 
ponent to the overall linewidth, simply by keeping track 
of how momentum conservation is achieved in each indi- 
vidual scattering event that contributes to the linewidth. 
In contrast, scattering of phonons or magnons from de- 
fects and boundaries is taken to contribute to the Umk- 
lapp component of the linewidth, as momentum is not 
conserved in such processes. Such scattering can often 
be extracted experimentally. 



In recent neutron spin-echo experiments on the 2D an- 
tiferromagnet Rb2MnF4, the combination of 4-magnon 
(two- in/two-out) scattering processes and boundary scat- 
tering was shown to agree with the magnon linewidth 
data[Tn] . Thanks to this excellent description, which was 
achieved without the use of adjustable parameters, nu- 
merical evaluation can be used to separate t„ and t„ for 
the 4-magnon scattering for the purpose of calculating 
Kniag within Callaway's framework. The boundary scat- 
tering was taken from the lowest-temperature data (at 
3 K)[Tn]. The main panel of Fig. 1 shows the resulting 
Kmag calculated for Rb2MnF4. This represents the first 
application of Callaway's theory without the use of any 
adjustable parameters. 

We also calculated the mean free path determined by 
the 4-magnon (two-in/two-out) interactions and the dis- 
persion relation (Fig. 1, inset). Below ^ 12 K, the 
mean free path is longer than the mean distance trav- 
eled to an intracrystalline boundary or defect, and the 
thermal conduction is thus determined by ballistically- 
transported spin-waves that collide with boundaries or 
defects, where they are absorbed and emitted. At higher 
temperatures, the conductivity is increasingly limited by 
magnon-magnon interactions within the magnon gas, and 
the boundaries play a less significant role. 

A qualitatively similar broad peak is also observed in 
the in-plane magnon-mediated thermal conductivity in- 
ferred from measurements in several square-lattice an- 
tiferromagnetic cuprates[SJ HJ |H]- We used the above 
method to calculate the thermal conductivity in the 
cuprate La2Cu04, which is isostructural to Rb2MnF4. 
Using the numerical methods described in Ref. 10[10). 
we calculated the linewidth for 4-magnon (two-in/two- 
out) scattering with 5 = 1/2 over the full Brillouin zone, 
using the low-temperature value for the nearest-neighbor 
exchange constant J in La2Cu04. We assumed that the 
renormalization of the exchange constant with tempera- 
ture in La2Cu04 follows that measured in Rb2MnF4 (ad- 
justed appropriately by S and 3 to reflect the different 
temperature scales). Because the magnetic anisotropy is 
considerably smaller relative to the exchange constant in 
La2Cu04 than in Rb2MnF4, for calculational purposes 
we treated La2Cu04 as if its dispersion were gapless; the 
calculated linewidth results thus represent a slight over- 
estimate. In these numerical evaluations, iterations were 
performed until the linewidths converged [lOj; this con- 
vergence was 0.6% or better over the full Brillouin zone 
for temperatures between 270 and 680 K. We took the 
boundary spacing L to be the temperature-independent 
mean free path of 558 ± 140 A that was fitted between 
70 and 158 K in Ref. 5 [5] and used this as the single 
adjustable parameter in calculating the component of 
the linewidth due to boundary scattering, as described 
in Ref. 10 [10]. 

In Fig. 2, we compare these calculated results for 
Kmag in La2Cu04 with representative extracted exper- 
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FIG. 2. Extracted and calculated values of Kmag in La2Cu04. 
The extracted values were determined by subtracting the esti- 
mated phonon contribution from the measured experimental 
data for k. 



imental results^ [H HH]. The agreement is reasonable, 
considering the strong dependence of the calculated re- 
sults on the boundary scattering (which was not mea- 
sured directly), and in light of the large variation be- 
tween samples IB HH]- Deviations in La2Cu04 from 
a simple Heisenberg model, including ring exchange T9J, 
may also be relevant. Furthermore, the calculated peak 
in La2Cu04 falls in the vicinity of Tjv, which is roughly 
325 K. At such high temperatures, the assumptions un- 
derlying the theory begin to break down, and other effects 
may become important. Finally, there are experimental 
uncertainties in the separation of the magnon-mediated 
and phononic contributions to the thermal conductivity. 

High-resolution magnon linewidth data over the full 
Brillouin zone thus establish a detailed correspondence 
between microscopic interactions (the magnon decay pro- 
cesses) and a macroscopic transport property (the bulk 
thermal conductivity). In order to evaluate the ther- 
mal conductivity from the linewidths of the excitation 
in question, quantitative description of the linewidths by 
theory is essential, as the contributions of normal and 
Umklapp processes cannot be separated experimentally. 
In Rb2MnF4, such agreement with theory was found (us- 
ing no free parameters) , and this agreement makes it pos- 
sible to link the magnon spectral function to the trans- 
port. The same approach can also be employed in phonon 
systems, for which theoretical descriptions of the thermal 
conductivity have heretofore involved multiple adjustable 
parameters 

This formalism provides a complete description of 
the thermal conductivity over a very wide temperature 
range. If the magnon (or phonon) dispersion relations 



are known, this description involves only a single un- 
known parameter, the domain boundary (or dislocation) 
distance L. This quantity can thus be extracted reliably 
from existing thermal conductivity data. This should 
provide a greatly improved basis for efforts to use the 
thermal conductivity as a probe of correlated metals and 
superconductors [H [2] and to design nanostructured ma- 
terials with thermal properties optimized for energy con- 
version and other applications [3]. 

The work in Stuttgart was supported by the German 
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